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ABSTRACT 


MODELING OF STRUCTURAL-ACOUSTIC INTERACTION USING 
COUPLED FE/BE METHOD AND CONTROL OF INTERIOR 
ACOUSTIC PRESSURE USING PIEZOELECTRIC ACTUATORS 

Yucheng Shi 

Department of Aerospace Engineering, Old Dominion University 
Advisor: Dr. Chuh Mei 

A coupled finite element (FE) and boundary element (BE) approach is presented 
to model full coupled structural/acoustic/piezoelectric systems. The dual reciprocity 
boundary element method is used so that the natural frequencies and mode shapes of the 
coupled system can be obtained, and to extend this approach to time dependent problems. 
The boundary element method is applied to interior acoustic domains, and the results are 
very accurate when compared with limited exact solutions. Structural — acoustic problems 
are then analyzed with the coupled finite element/boundary element method, where the 
finite element method models the structural domain and the boundary element method 
models the acoustic domain. Results for a system consisting of an isotropic panel and 
a cubic cavity are in good agreement with exact solutions and experiment data. The 
response of a composite panel backed cavity is then obtained. The results show that the 
mass and stiffness of piezoelectric layers have to be considered. 

The coupled finite element and boundary element equations are transformed into 
modal coordinates, which is more convenient for transient excitation. Several transient 
problems are solved based on this formulation. 

Two control designs, a linear quadratic regulator (LQR) and a feedforward controller, 
are applied to reduce the acoustic pressure inside the cavity based on the equations in 
modal coordinates. The results indicate that both controllers can reduce the interior 
acoustic pressure and the plate deflection. 


FOREWORD 


The research results contained in this technical report were performed under a NASA 
Langley Research Center Grant NAGI-I684 entitled, "Experimental and Numerical Analysis 
of Structural Acoustic Control for Interior Noise Reduction." The report is based essentially 
on a Ph.D. dissertation prepared by Yucheng Shi under the supervision of Dr. Chuh Mei. The 
document presents a coupled finite element and boundary element formulation for the modeling of 
fully coupled structural/acoustic/piezoelectric systems and the active control of interior acoustic 
pressure using piezoelectric actuators. The work was conducted at the Department of Aerospace 
Engineering, Old Dominion University. Mr. Travis L. Turner of Structural Acoustic Branch, 
NASA Langley Research Center, was the technical monitor. 
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Chapter 1 
INTRODUCTION 


1.1 Preliminary Remarks 

The physical manifestation of sound is a time-dependent pressure fluctuation around 
the static pressure in a compressible fluid, such as air or water. One such source for 
these pressure fluctuation is a vibrating elastic structure. The fluctuating pressure on 
the surface of the structure constitutes the radiation loading. Generally, because of 
the low density of air compared to structural materials, radiation loading exerted by the 
atmosphere is usually small enough to have a negligible effect on the structural vibrations. 
Consequently, the two theoretically coupled systems, the elastic structure and its dynamic 
response within the atmosphere under prescribed driving forces and the acoustical pressure 
field generated by the velocity distribution over the structure-atmosphere interface, can 
be analyzed independently. 

However, this is not always the case. For example, when a volume of air in 
contact with the structure is confined in a small enclosure, or when the structure is 
exceptionally light, the influence upon the structural motion due to radiation loading has 
to be considered. In those circumstances, the structure vibrating in air is more like a 
structure vibrating in contact with a fluid of comparable density, where radiation loading 
is comparable to the inertial and elastic forces of the structure. Hence, the elastic and 
acoustical dynamics, and their interactions, must be modeled simultaneously. 

True structural-acoustic systems can be found in numerous industrial applications 
such as interior noise, or noise transmitted into a cavity. For example, the sound inside 
a fuselage of an aircraft with engines or fluid flow as outside sound sources, or the noise 
inside the automobile passenger compartment with the engine or driveline vibration or 
body vibration as the outside sources are important engineering problems. Shown in 
Figure 1.1 is a general sketch of these types of problems. The sound sources and the 
cavity are separated by an elastic structure. The vibration of the elastic structure, here a 
flexible plate, is excited by radiation loading on the exterior surface of the plate due to 
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the sound source. This induces a sound pressure fluctuation field inside the cavity due 
to the velocity distribution on the inner surface of the plate. The coupling between the 
acoustic field and structure vibration on the outer surface is negligible because of the low 
density of air compared to the structural material. Radiation loading on the outer surface 
of the flexible plate acts like an external force only and the plate has no effect on the 
outside acoustic field except reflection and absorbing. On the other hand, the coupling 
between the plate vibration and the acoustic field inside the cavity has to be considered. 


Acoustic Source 



Figure 1.1 The structural-acoustic interaction problem 

Over the past several decades, aerospace and automotive industries have given more 
attention to acoustic excitation, because this is not only a passenger comfort issue but 
also a safety concern. The acoustic pressure level outside an aircraft fuselage can reach 
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164 dB [11 , and may go as high as 190 dB [2] for some new supersonic aircraft, indicating 
radiation loading on the fuselage is an important design parameter. Under such conditions, 
the interaction between the internal acoustic field and the structure vibration becomes 
extremely important. Mathematical models which can predict the response of the coupled 
structural acoustic system under prescribed external forces becomes a necessary tool 
for designers. 

The prediction of the sound pressure field inside the cavity and the dynamic response 
of the plate is not the final goal of analysts. Other objectives are to reduce the inside 
pressure level and the structure vibration. These objectives are not independent of the 
modeling of the structural-acoustic interaction. In addition to the requirement for a control 
system, more features are required in the coupled structural — acoustic mathematical 
model to represent the control design based on piezoelectric materials. 

1.2 Review of Previous Work 

Research in both the modeling of the structural-acoustic interaction problem and the 
active control of structurally transmitted noise has been performed. In this section, a 
brief review of both topics is given. 

Modeling of Structural-Acoustic Interaction 

It is not until recent years that the structural-acoustic interaction problem has captured 
interest in the aerospace and automotive industries. However, the basic problem is an 
old one. The oldest studies can be traced back to World War I, when Rayleigh published 
the first modem text on acoustics, the Theory of Sound^. He formulated the equation of 
motion of a rigid spring piston radiating into an acoustic fluid, and considered the effect 
of the acoustic fluid by increasing the damping and mass of the single-degree-of-freedom 
system.^ 1 

Even the study of radiation of sound by a vibrating structure into an acoustic cavity 
has a history of more than 30 years. The first investigation was performed by Lyon [5] in 
1963. In his work, a rectangular plate backed by a rectangular cavity was studied in a 
straightforward but approximate manner. This problem was then investigated by many 
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researchers such as Dowell and Voss 161 , Bhattacharya and Crocker 171 , Guy 18, 91 , McDonald 
et al. 1 In 1977, Dowell, Gorman, and Smith 1 ' ' 1 derived the general governing equation 
for the coupled structural-acoustic system, so-called “acoustoelasticity”. In this classic 
paper, a comprehensive theoretical model was developed for interior acoustic fields 
which were created by flexible wall motion resulting from exterior sound fields, and 
accurate coupling between the wall and interior acoustic cavity was considered. A modal 
interaction approach, which assumed that the coupled modes of the cavity wall system can 
be expressed as a linear combination of the rigid wall cavity acoustic modes and invacuo 
plate modes, was also proposed for multiply connected cavities. The comparison of the 
numerical results and experiment data was shown to be in good agreement 11 Based on 
Dowell’s formulation, extensive works by others in the field followed 112-161 . Bokil 1171 
obtained a closed-form solution for the acoustic pressure transmitted through a rectangular 
cavity backed flexible plate. 

On the other hand, in 1966, Gladwell and Zimmermann 1181 developed an energy 
formulation of the acoustic-structure interaction problem, this paper set the stage for 
the application of finite element methods to cavity-structure analysis. This numerical 
method makes the consideration of complex cavity and structure geometry, structure 
boundary condition, and acoustic boundary condition conceptually no more difficult than 
simpler problems. Three different formulations were derived using the pressure 119-231 , 
fluid particle displacement 124-271 , or velocity potential 1281 as the fundamental unknowns 
in the fluid region. The finite element approach for the structural-acoustic interaction 
problem seems well developed. In 1970’s, even the computer tool NASTRAN had 
the capability of cavity-structure analysis 129,301 . Neffske et al. 1311 analyzed the acoustic 
pressure field of complex automobile passenger compartments using NASTRAN in 1980 
and found good agreement with experiment data. 

As a powerful alternative to the finite element method, the boundary element method 
(BEM) or the boundary integral element method (BIEM) had its beginnings in the early 
1960s based on the boundary integral equation theory developed in 1800s and 1900s 1321 . 
This method was first applied to the acoustic area to solve an acoustic radiation problems 
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by Chen and Schweikert in 1963 [33] , followed by extensive research work from 1960s to 
1990s [ 2 \ Most of the boundary element method applications in acoustics focused on 

the acoustic radiation and scattering problems, where boundary element methods have an 
incompatible advantage for dealing with infinite domain. Not until 1982 did Koopmann 
and Benner present the application of this method to internal domains 143 !. Suzuki et al. 144 !, 
Tanaka and Masuda! 45 - 46 !, Mariem and Hamdi [47] , and Pates 148 ! applied this method 
to structural-acoustic interaction problems assuming a sinusoidal time dependence. A 
common feature of those investigations is the joining of the boundary element method 
with the finite element method. The finite element method was used to model the 
structures, while fluid domain was handled by a boundary element method. 

The most important issue for boundary element methods is the selection of the 
fundamental solutions. For an internal acoustic problem, the general governing equation is 
the well-known three dimensional wave equation. Based on various different assumptions, 
different choices of the fundamental solution lead to different approaches 140 !. If the 
fundamental solution of the wave equation is selected 150 !, the boundary element equations 
need no domain integral, but it becomes impossible to form an eigenvalue problem to 
obtain eigenvectors. If the fundamental solution of the Laplace operator is used, the 
time-dependent term will inevitably lead to domain integrations 151 !. 

Active Control of Acoustic Pressure Level Inside the Cavity 

There are two different concepts in active noise control, one is the use of secondary 
sound sources, such as microphones, to reduce the undesired noise 152-63 !, the other is to 
reduce the original noise source 164-68 !. For transmitted acoustic pressure through elastic 
structure, using piezoelectric materials embedded in the structure to reduce the pressure 
level belongs to the later case. The piezoelectric actuators will control the vibration of 
the elastic structure so that the acoustic pressure level transmitted inside the cavity can 
be reduced. Note the objectives of active control may involve more than just reduction 
of structure deflection. 
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Early attempts at reducing radiated sound from a vibrating structure by active control 
utilized point forces in a feedforward control loop [57,58] . This approach was then extended 
to applying bending moments and in-plane forces using piezoelectric actuators 164 ’ 661 . 
Analytical and experimental work supported the feasibility of using active vibration 
control through force inputs in structures to reduce interior cavity noise level 165, 681 . 
However, most all studies have been based on modal analysis or frequency response 
input/output analysis, which does not include direct coupling of the structure and the 
acoustic cavity 1681 . Banks et al . 166, 671 proposed a time-domain state space formulation 
based on the finite element model of a two-dimensional rectangular cavity backed by 
an isotropic beam. The structural-acoustic interaction was considered, but the coupling 
between the structure and the piezoelectric actuators was not complete. Only the control 
force from the actuators was included, the mass and stiffness of the actuator, or the 
piezoceramic patches, were ignored. 

On the other hand, the control system design based on the coupled structural- 
acoustic model was performed using secondary sound sources 152,541 . Snyder et al . 1561 
presented a theoretical framework suitable for control system design. They proposed 
various control objectives such as minimization of acoustic potential energy in the cavity, 
minimization of acoustic pressure amplitude at a discrete location, or minimization of 
structure kinetic energy. Thomas et al. 160, 611 presented a numerical control simulation of 
sound transmission through a cylindrical shell using secondary sources. Results indicated 
the possibility of reducing the acoustic potential energy in the cylinder and the acoustic 
pressure amplitude at a discrete location using various weighting matrices in a linear 
quadratic regulator controller. 


1.3 Objective and Outline 

The overall objective of the present study is to apply available control techniques 
to reduce the transmitted acoustic pressure level inside the cavity. To perform this task, 
three steps are required. The first objective is to develop a coupled finite element and 
boundary element method to model the structural-acoustic system. The second step is 
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to transform the coupled equations into the time-domain modal formulation. The third 
step is to apply the linear quadratic regulator and least-mean-square (LMS) controllers to 
reduce the sound pressure level inside the cavity based on the modal formulation. 

Since the control design is based on the numerical model, some restriction must 
be considered in the structural-acoustic model. The finite element method used to 
model the flexible plate has to be able to deal with non-homogeneous materials, that 
is, a composite laminate with embedded piezoelectric material layers. This involves 
not only the mechanical and electric coupling properties of piezoelectric materials, but 
also the mass and stiffness of the piezoelectric layers which have never been considered 
in studies reported in the literature. This feature can be easily accommodated with 
the versatility of the finite element method, but the requirement of transient response 
predictions on the combined finite element and boundary element method for the coupled 
structural-acoustic system does need some additional effort, specially for the boundary 
element model of the acoustic cavity. The boundary element method used has to be 
able to handle the time-dependent loads and also be able to calculate the eigenmodes 
of the coupled system. The dual reciprocity boundary element method^ 69, 7 °J is used in 
present study, this new coupled finite element and boundary element method introduces 
an innovative and powerful approach. The formulation of the finite element method 
for a composite plate with embedded piezoelectric material layers, the dual reciprocity 
boundary element method for the time-dependent acoustic problem, and the coupled finite 
element and boundary element method for the structural-acoustic interaction problem are 
all developed in Chapter 2. The numerical results and their comparison with exact 
solutions, experimental data, and other numerical approaches are given in Chapter 3. 

The modal formulation is also given in Chapter 2 with numerical results for various 
loading cases presented in Chapter 3. To transform the coupled equation into modal 
formulation, the modal coordinates, or the basis of the solution space, has to be determined 
first. The modal shapes of the coupled structural-acoustic system are used here rather 
than the uncoupled modal shapes. One of the advantages of this treatment is that the 
basis will yield a set of uncoupled differential equations, instead of coupled equations. 
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The control portion of this study is to demonstrate that control laws can be applied to 
reduce acoustic transmission effectively. The linear quadratic regulator shows how much 
one can reduce the transmitted sound level in a theoretical sense. The feedforward least- 
mean-square control law is an attempt to design a more practical controller. The control 
formulations and numerical results are all given in Chapter 4. Concluding remarks and 
recommendations for future work are presented in Chapter 5. 
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Chapter 2 

COUPLED FINITE ELEMENT AND BOUNDARY ELEMENT 
METHOD FOR STRUCTURAL-ACOUSTIC INTERACTION 

In this chapter, the finite element method for a composite plate with embedded 
piezoelectric layers is developed, and the dual reciprocity boundary element formulation 
for acoustic problem is derived. Following this treatment, the finite element method and 
boundary element method are coupled together to form the discrete governing equations 
for the structural-acoustic interaction problem. 

2.1 Coupled Finite Element Equation of 
Motion of Composite/Piezoelectric Materials 

The finite element method has become a very powerful tool in the analysis of static 
and dynamic response of structures. This method is capable of handling structures having 
complex geometries, non-uniform materials, and complex boundary conditions. The 
objective here is to develop a finite element model to describe the coupling between a 
plate or structure and piezoelectric materials. 

2.1.1 Piezoelectricity 

It is well known that piezoelectric materials can be used as actuators and sensors 
due to their direct and converse effects 17 ’ 1 . To model these effects, a formulation 
with mechanical and electrical coupling should be considered. This electromechanical 
formulation presented here is based on the linear piezoelectricity theory 172,731 . In this 
section, the general piezoelectric constitutive equations are briefly reviewed, and related 
piezoelectric constants are defined. 

The electric enthalpy H is defined as the amount of energy stored in the material 1721 

H = U — E • D (2.1) 

where U is the total internal mechanical energy, E and D are the electric field and 
electric displacement vectors, respectively. Let a and e be the mechanical stress and 
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strain tensors, they are related through the electric enthalpy as [72] 


_ dH 

a ~ de 

Similarly, we have the relationship between the electric field 
vector as^ 72 ' 


( 2 . 2 ) 

vector and displacement 


D 


dH 

dE 


(2.3) 


The electric enthalpy H is assumed to have a homogeneous quadratic form in the 
linear piezoelectric theory, which includes the contributions from elastic strain energy, 
piezoelectric energy, and electric energy and is given by 

. t _ . . 

(2.4) 


^ r,) rmTr^ir.i 1 r m T r 


H = §W' MW - { E}‘ [e]{e} - j{ E}‘ [=]{£} 


where superscript T represents transpose, [<J] is the stiffness matrix assumed independent 
on electric field, [c?] is the dielectric matrix assumed independent on strain, and [e] is 
the piezoelectric constant matrix which relates stress to applied electric field, and that is 
where the coupling electrical and mechanical features originate. 

Substituting equation (2.4) into equations (2.2) and (2.3), the linear constitutive 
relations become 


M = MW - M t {£} (2.5) 

and 

<D}=[e]< t } + M r {£} (2.6) 


Since the dielectric matrix [w] and the piezoelectric constant matrix [e] are sometimes 
not available, the more commonly available constant matrix [d] and the free dielectric 
matrix [e] are introduced, the relationship between those matrices are^ 73 ^ 


[e] = [d\[Q] (2.7) 

and 

[e] = [w] + [d\[QM T (2.8) 

The constitutive equations of piezoelectric materials become 

M = W](M-M t {£}) (2-9) 

(£)) = [<(lH + [£] 7 '{£) (2.10) 
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2.1.2 Constitutive Equations for Composite and Piezoelectric Laminae 

For a thin laminated panel consisting of composite and piezoelectric layers (shown 
in Figure 2.1), the two-dimensional constitutive equation in principle material axes for 
an orthotopic piezoelectric layer can be obtained from equation ( 2 . 9 ) as [71] 
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( 2 . 11 ) 


and the electrical displacement along the polling 3 -axis, which is assumed to be the 
normal direction of the plate, from equation ( 2 . 10 ), is 


D$ — [cfoi d$2 0] < 



( 2 . 12 ) 


where is the electric field; d^i , d ^2 are the piezoelectric stress/charge constants; £33 is 
the permittivity constant. For a composite lamina, the constitutive equation is 
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( 2 . 13 ) 


where [Q] is the reduced stiffness matrix, and the subscripts <j> and c denote for piezo- 
electric and composite laminae, respectively. 
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Figure 2.1 Composite and Piezoelectric Layers 
Then the strain-stress relation for a general k-th layer with a lamination angle becomes 
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where for composite layers, 

[ 0 ] = [$]. 

E 3 k = 0 

and for piezoelectric layers, 

[ 0 ] = [< 31 * 

E 3 k = (d x , dy, [Q] ^ {e} + E 3k £ 33k 

where [0] is the transformed reduced stiffness matrix. 


(2.14) 


(2.15) 


(2.16) 
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2.13 Finite Element Equations 

The generalized Hamilton’s principle for a coupled electromechanical system can be 
expressed as^ 74, 75 1 


l 2 

j (T-U s + W e -W m + W)dt 


Ui 


= 0 


(2.17) 


where 


(2.18) 


(2.19) 


T 1 f dwdw 

t= 2 J p -mm dv 

V 

is the kinetic energy; 

Us = \j W T ^ dV 

v 

is the strain energy of the system; 

W' = \ J {E) T {D)dV 

V 

is the electrical energy; W m is the magnetic energy, which is negligible here; and 

W = J wF b dV + J wF s dS + wF c - j <f>qdS (2.21) 


( 2 . 20 ) 




is the work done due to external forces and the applied surface charge. In equations 
(2.18) to (2.21), p is the density, w the displacement vector, F b denotes the body force 
vector, F 9 the surface traction vector, F c the concentrated forces, <fi the electric potential, 
and q the surface charge. 

The linear strain-displacement relation is given as 


( 2 . 22 ) 
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and the relation between electric field and electrical potentials is 

{£ 3 } = -[B h \{4>) 


(2.23) 
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where [F^] is a diagonal matrix 


1 -l 


[B h } = 


'h i 
0 



(2.24) 


with h k is the thickness of k-th piezoelectric layer, u, v are the in-plane displacements, 
w the transverse displacement, and is the number of piezoelectric layers. 

Substituting equations (2.18) through (2.21) into equation (2.17), and using equations 
(2.14) to (2.16) and (2.23) to (2.25), we have the equation of motion for the composite 
plate with embedded piezoelectric layers in finite element form as 


ch::} 

where {W} and {$} are the system structural node degrees-of-freedom and electric po- 
tentials, {F w } and { F^ } are the load vectors due to acoustic excitations and piezoceramic 
surface charge, and [A/] and [I<] are the system mass and stiffness matrices. The sub- 
scripts w and <f) denote structural and electric field components, respectively. The detail 
derivation, element load vectors, and mass and stiffness matrices are given in Appendix A. 
The equations of motion can then be rewritten as actuator equations as 

{Mjlw} + \K w \{W} = {F„,} + 

(2.26) 

{fW} = -[iC*]W 

and as sensor equations as 


M w 0 
0 0 


VV 


+ 


K w 

yK<j>w 


K w4 > 

K<t> 


[K*] {^} + [K^] {W} = {F^} (2.27) 

where {F^} is the load on structure due to electric potential {$}. 

2.1.4 Finite Element Specifications 

Many rectangular and triangular type finite elements are currently being used in 
commercial and in-house codes. Any type of finite element can be applied to the present 
formulation. The element selected for this study is a four-node rectangular C 1 conforming 
element 176 ! anc j } s shown in Figure 2.2. Each node of this element has four degrees of 
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freedom associated with the transverse displacements (w,w jX ,w t y ,w <xy ) and two in-plane 
displacements (u, u). The electric degrees of freedom are the electric potential <j> of each 
piezoelectric layer. The displacement shape functions and other characteristics of this 
element are given in Appendix B. 



Figure 2.2 C 1 conforming rectangular finite element 


2.2 Boundary Element Method for Acoustics 


2.2.1 General Aeroacoustics Equations 


The general governing equations of inviscid flow are formed by one continuity 


equation [77 ^ 




0 


(2.28) 


three momentum equations 


and an energy equation 



dp - , 

+ dZ =iifi 


(2.29) 
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where c is the sound speed, p is the density of the medium, u{(xj,t), i,j = 1,2,3 are 
the velocity components, p(xj,t) is the pressure, and /, are body forces. 


The change of density and pressure due to acoustic wave is usually much smaller 
than that due to fluid flow, so it is necessary to separate acoustic wave from the fluid 
flow, denote subscript ’o’ as the variable of fluid flow (i.e., p = p + p 0 ), the perturbation 
decomposition gives the governing equations of acoustic wave in a homogeneous medium 
with no mean flow as 


dp dui 

dui dp 

^ + d 7 t =° l = 1 ’ 2 ’ 3 

^-c 2 ^ = 0 
dt dt 


( 2 . 31 ) 


Eliminating the perturbation of density from above equations, we have the governing 
equation as 


d 2 p 

dt 2 


- c 2 v 2 p = 0 


( 2 . 32 ) 


and the momentum equations give the relationship between acoustic pressure and velocity 


components as 


dp du{ 

§^ + ' , ' > ar = 0 


( 2 . 33 ) 


and at a boundary, the normal derivative of pressure is obtained as 


9 


dp duiUi 

dn dt 


( 2 . 34 ) 


where n is the unit outward normal vector. 

Generally, it is more convenient to use a velocity potential to simplify the equations, 
but here, since both pressure and velocity are needed in the formulation of control theory, 
we use the pressure as the primary variable. 


2.2.2 Dual Reciprocity Method 

As one type of the boundary element method, the dual reciprocity method^ 69, 70 ^ 
offers the advantage that it does not need any domain integrals which appear when 
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nonlinear problems or time dependent problems are considered (which does not mean 
no domain equations are needed). However, in this study, the main reason for using 
the dual reciprocity method instead of the general boundary element method is that we 
have to transform the coupled structural — acoustic equations into modal coordinates. To 
do so, some kind of uncoupled coordinates are required, and one commonly used set of 
coordinates are the eigenvectors. Generally, the boundary element method combines 
the frequency of the input excitation into its fundamental solution 1481 , making the 
identification of eigenvalues and eigenvectors extremely difficulty. 


General formulation For a general boundary value problem, 

V 2 p(x,y,z,t) = b(x,y,z,t) in ft (2.35) 

Assume that the right-hand term, b, can be expressed as a linear combination of a set of 
independent functions, /,, defined in the domain ft and on the boundary T. That is, 

NN 

& = £>,/, (2.36) 

i=i 

where a, are initially unknown coefficients, and NN is the number of independent 
functions. Defining a set of particular solutions, pi, as 

V 2 p, = /, i = 1, 2, ..., NN (2.37) 

Then equation (2.35) becomes 

NN 

V 2 p=£>,(V 2 p,) (2.38) 

t=l 

If we take the fundamental solution of the Laplace’s equation 

VV(f) + <5(r J ) = 0 (2.39) 


that is 
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AkR' 
dp * 

dn 


R= \r-fj\ 

1 dR 
4 nR 2 dn 


(2.40) 
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as the weighting function, then the weighted-residual method applied to equation (2.38) 
requires 


J (v 2 p)p*dn = J2<** f (v 2 Pi) P *dn 

O 1 — 1 o 


(2.41) 


n 1=1 n 

Using Green’s formula, we obtain an integral equation as 


J(V 2 p’)pdSl + J (pq'-p-q)<tr = 

« r 



where T is the boundary of the domain. 


Using equation (2.39), and considering the singularity of the fundamental solution at 
the boundary, we obtain the integral equation for each point on the boundary as 

{m* ~ (2.43) 

where parameter c } is the ratio of the exterior area outside the domain of a sphere surface 
with a small radius to the area of the entire sphere surface. For any interior point, we 
have cj = 1 and 



( m * ~ P*4i)dY J (2.44) 

Equation (2.43) is the basis of the boundary element method. Since we seek a time- 
dependent formulation, equations should be formed for nodes in the domain and on the 
boundary, but as can be seen from equation (2.44), no domain integrations are required. 
If the boundary is discretized into elements, and in each element, shape functions are 
assumed so that the variables defined on the boundary can be expressed in terms of their 
node values or 



P = [WHp} 


then the integral equation becomes 


N N NN T 


- ^2 HjkPk - ^2 GiMk = *22 

jt=i jt=i i=i 


<*« 


N N 

£ - E 

Jt=l Jt=l 


(2.45) 


(2.46) 
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where [iV] is the shape functions, and N is the number of boundary elements. In matrix 
form, we have 


-WW-IGIW) = (-(#][£] -[G] [«])w (2.47) 


where matrices 


and 


Q 


contain p, and q, as their columns, respectively, matrices 
[#] and [G] are globe acoustic influence matrices assembled from the corresponding 
element matrices which are defined as 


[sli = J p‘\N]<ir 

r* 

W k = -cjS }k - J g *[yv]dr 

r fc 


(2.48) 


Functions As expressed in equation (2.36), when function b satisfies certain continuous 
conditions, there will be a set of functions and correspond coefficients such that 


{6} = [F]{a] (2.49) 

where each column of [F] consists of a vector which contains the value of the function 
f, at the (NN) collocation points. Then if matrix [F] is invertible, we have 

M = [^T 1 ^} (2.50) 


Substitute the above equation into equation (2.47), and denote 


we obtain 


[S] = (-[tf] 



!G][o])(fr‘ 


(2.51) 


-|»l(f] - [G]{<?} = [SUM (2.52) 

Matrix [5] is a function of /,, and p, only, where the requirement on functions /„ 
except some continuous requirements, is to allow that matrix [F] to be invertible. One 
choice for those functions is to let ^ 78, 79 ^ 


/.(f) = 1 + R„ R, = |f - f,-| 


(2.53) 
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and for a three-dimensional problem, we have 


„ R* R} 

Pi = — ~ -| 

12 6 


(2.54) 


and 


s , Ri'SdRi 

q, ~ an “VT + T)as 


(2.55) 


When these three sets of functions are obtained, matrices [H], [G] and [S] are known, 

and with the function b, equation (2.52) can be solved with suitable boundary conditions 
and initial conditions. 


2.23 Boundary Element Method Equations 


For a time-dependent acoustic problem, equation (2.32) is the governing equation, 
comparing to equation (2.35), we set 


b = 



(2.56) 


so from equation (2.52), we have 


m{P} +[/(]{/>} + [(?]{(?} = 0 (2.57) 


where 


Wa) = ^[5] (2.58) 

2.2.4 Boundary Element Specifications 

Similar to the finite element method, many types of boundary elements are available. 
In this study, we selected one of the simplest two dimensional elements, the constant 
rectangular element (shown in Figure 2.3), to demonstrate the capability of boundary 
element method in solving acoustic problems. 
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Node Representing Element 



Figure 2.3 Two dimensional constant boundary element 


2.3 Coupled Structure/Acoustic Equations 

The boundary element equations of acoustic wave inside the cavity can be rewritten 
here as 


M c M c b 

Mbc Mb 



H c H cb 
H hc H b 



, f °c G cb 

G bc G h 



where the pressure and its normal derivative are divided into two parts. 



(2.59) 


(2.60) 


the subscript c denotes the quantities on the inner surface of the structure which will be 
coupled with the finite element structure equation, and b at all the other locations. The 
governing finite element equation of the structure, equation (2.26), is 


[M w }{w] + [K w ]{W) = {F ei } - {F in } + {F^} (2.61) 
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where the acoustic load is separated into external excitation load and inside acoustic 
load, {F w } = {F ex } - {F ln }. 

Those two set of equations will be coupled by using the compatibility of displacement 
and pressure on the inner surface of the plate. 

2.3.1 Coupled Equation in Structure DOF 

Because of the continuity of the displacement, the normal component of acoustic 
velocity on the inner surface of the plate must be equal to that of the plate at the same 
location, or 


{w„c} = [^]{uc} = [/?][T]{^} (2.62) 

where the matrix [/?] transforms the acoustic velocity on the inner surface of the plate 
{t/ c } to its normal component {u nc }, and matrix [F] is a transformation, which depends 
on the primary variables and meshes of finite element and boundary element models. The 
element matrices [/? e ] and [F e ] for constant boundary elements are derived in Appendix 

C. Since equation (2.34) gives a relation between pressure and acceleration components, 
we have 


{Qc} = -Po{<c) = -Po[P][T}{w} = [7i]{if} (2.63) 

The second loading term on right side of equation (2.61) is due to the acoustic 
pressure inside the cavity, its element node value can be related to the inside acoustic 
pressure (see equation A. 16) as 

{fin} = [T b ] T f Pc {H w }dA (2.64) 

A 

Form this, we can obtain the relationship on a global level as 

{fi„} = Y. (inf I {«A<ia \ {/>,} = [r 2 ){P c ) (2.65) 
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Substituting equations (2.63) and (2.65) into equations (2.59) and (2.61), we have 


the coupled equation for a structural — acoustic system as 
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[M](x} + [C]{x} + [AIM = {Ai} + {F 2 } 
where {x} = [W,P c ,P b ] T and 

([M cw ),[Mbw]) = ([Gc},[G b c})[Ti} 

[K wc ] = [T 2 ] 


( 2 . 66 ) 


(2.67) 


( 2 . 68 ) 


Notice that the generalized mass matrix [M] and the generalized stiffness matrix [A] 
are all no longer symmetric, and for generality, a damping term [ C w \ j VK j is included 
in equation (2.66). The force {Fi } consists of external excitation to the structure and 
normal derivative of pressure in the cavity, and {A 2 } is the force acting on structure due 
to applied electric potential. 


2.3.2 Coupled Equations in Modal Coordinates 


There are at least two different modal coordinate systems that can be used to transform 
the governing equation in node DOF to truncated modal coordinates. One is taking the 
natural modal shapes of the uncoupled system, defined by equation (2.59) and equation 
(2.61). This allows the nodal DOF to be expressed as 


{x} = 





(2.69) 


where ['M is the modal matrix for plate equation, ['£ a ] ([^a] = [^c, "f&jy is the 
modal matrix for acoustic equation, and { 77 ^} and {ri a } are the modal variables of plate 
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displacement and acoustic pressure, respectively. Another transformation is based on the 
truncated natural modal shapes of the coupled equation (2.66), which gives 

{*} = Wiv) (2.70) 

where { 77 } is the vector of modal variables, and [\D] is the modal matrix. 

The first transformation will end with a set of coupled ordinary differential equations 
because of the coupling between structure and acoustics even when no damping is 
involved. Further, cases with all rigid acoustic boundary condition require some special 
treatment for the rigid acoustic mode. Meanwhile, although more computing time is 
needed to form coupled modes, the second transformation will give a set of independent 
ordinary differential equations in the case without damping, and no judgement is needed 
to select modes. Therefore, the second transformation, using coupled modes as modal 
coordinates, is applied here. 

Substitute equation (2.70) into equation (2.67), and pre-multiply by the transpose of 
corresponding modal matrix, or 

+ [crw + m = (m + { l 2 } 

where the modal matrices and modal vectors are 

««]MCT,l*n = (¥] r ([Afl,(C], [*•])[*] 

({M.{£2}) = M r ({Fi},{f2}) 

Equation (2.71) is the so-called time-domain modal formulation. 


(2.71) 


(2.72) 
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CHAPTER 3 

RESULTS FOR THE MODELING OF 
COUPLED STRUCTURAL-ACOUSTIC SYSTEM 


In this chapter we present the results based on the formulation previously developed 
and also some additional discussion. The results are subdivided into three sections: acous- 
tic results using the dual reciprocity boundary element method, the coupled structural- 
acoustic response using the coupled boundary element and finite element method, and 
results using the time-domain modal formulation. The purpose of this chapter is to show 
that the dual reciprocity boundary element method is a powerful numerical technique that 
can be used to solve a wide range of acoustic problems. The coupling of the boundary 
element method and the finite element method, and its time-domain modal formulation, 
lead to powerful tools that can be readily used to analyze structural-acoustic interaction 
problems. 


3.1 Interior Domain Analysis 

In this section, three example problems are solved and results are compared with 
exact solutions. Those problems are derived from a typical rectangular duct as shown in 
Figure 3.1. The inlet is assumed to be at 2 = 20cm and the outlet at z = 0, while the 
four duct side walls, with length of 20cm, are all assumed to be acoustically rigid in all 
three examples, that is = 0. The outlet has different boundary conditions. 

The input acoustic wave, that is the acoustic field at the inlet, is assumed to be a 
unit plane wave, thus the one-dimensional exact solution becomes valid 1801 . While for 
the dual reciprocity boundary element approach, the inlet, outlet, and the four duct walls 
are all discretized by constant two — dimensional boundary elements, and internal nodes 
are uniformly distributed inside the duct. Two different models with different mesh are 
employed for each problems. The first model (model- 1) discretizes the inlet, outlet, or 
each wall by a 5 x 5 mesh (or 150 boundary elements for the whole boundary) and 
5 x 5 x 5 (or 125) internal nodes. The second model (mode-2) discretizes the inlet, 
outlet, or each wall by a 7 x 7 mesh (or 294 boundary elements for the whole boundary) 
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and also 5 x 5 x 5 (or 125) internal nodes. In the following example problems, the 
pressure obtained at location x = y = z = 10 cm is compared with the exact solution 
at different input frequencies. 

In the first example, a null pressure field is applied at the duct outlet, the results are 
given in Figure 3.2. It can be seen that the boundary element method gives very accurate 
results compared to the exact solution. However, at the acoustic natural frequency of the 
duct, the boundary element results vary from the exact solution slightly. The source of 
this discrepant is the coefficient matrix singularity at this frequency. 

The second example problem uses the same duct and the same acoustic assumption, 
but the outlet boundary condition is assumed to be a rigid wall. Once again, from 
Figure 3.3, the boundary element method gives very accurate approximation, except at 
the acoustic natural frequency. 

In the third example problem, an impedance boundary condition is applied at the 
outlet of the duct. A non-reflection condition, that is no reflection from the outlet, is 
used here, which assumes the relationship between pressure and velocity at the outlet 
to be v = p 0 c based on the plane wave assumption. The density of air is taken as 
p 0 — 1.21 kg/m and the speed of sound in air is given as c = 343 m/sec. The results 
are shown in Figure 3.4. Compared to the exact solution^ 481 , observe that the boundary 
element method results are quite accurate. 

Generally, when the input frequency increases, the boundary element method results 
tend to diverge from the exact solution. At high frequency, the higher acoustic modes 
become significant while the current boundary element mesh may not be fine enough to 
describe the higher acoustic modes. Thus, a finer mesh is needed to resolve the solution. 


26 





Frequency (Hz) 


100 9< 


Figure 3.2 Results for three-dimensional duct with null pressure at outlet 








Pressure 



Figure 3.4 Result for three-dimensional duck with non-reflection outlet 
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3.2 Structural — Acoustic Interaction Response by Coupled FE/BE Method 

Coupling the structural and acoustic domains is essential when trying to accurately 
model harsh environments. As mentioned in Chapter 1, the coupling between a flexible 
plate and the acoustic cavity under harmonic acoustic inputs has been investigated by 
many researchers using different approaches. The first objective of this section is to 
validate the present coupled FE/BE method by comparing with known analytical and 
experimental results for a system consisting of an isotropic plate and a cubic acoustic 
cavity. Then this coupled FE/BE method is extended to other systems with composite 
plates. 

As shown in Figure 3.5, the coupled structural-acoustic system studied in this section 
consists of a square plate and a cubic acoustic cavity with all other five walls acoustically 
rigid. The plate is subjected to a uniformly distributed external pressure loading. Two 
different plate models are considered. 

The first coupled structural — acoustic system studied here consists of a cubic cavity 
and a simply supported brass panel with the following plate and cavity characteristics: 

Panel/Cavity Problem-Simply Supported Brass Plate 

Cavity 

X axis length L x = 20 cm 

Y axis length L y = 20 cm 

Z axis length L z = 20 cm 

Density of Air p 0 = 1.21 x 10~ 3 g/cm 3 
Sound speed c = 34300 cm /sec. 

Brass 

Young’s Modulus = 10.4 x 10 10 psi 
Poisson’s Ratio = 0.37 
Density pb = 8.5 g/cm 3 
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Brass Plate 


X axis length = 20 cm. 

Y axis length = 20 cm. 

Thickness = 0.09144 cm. 

Two different mesh cases were used in modeling this coupled system. One case 
(mesh-1) consisted of 25 (5x5 mesh) finite elements for the plate, 150 (5x5 mesh for 
each of the six cavity walls) boundary elements for cavity, and 125 (5x5x5) internal 
nodes; the second case (mesh-2) used 49 (7x7 mesh) finite elements for the plate, 294 
(7x7 mesh for each of the six cavity walls) boundary elements, and 343 (7x7x7) internal 
nodes. The pressure at the center of the cavity back wall ( at x = y = 10 cm, 2 = 0) is 
considered as the output pressure (p ou t ) and the transmission loss (TL) is calculated using 
the pressure on the exterior surface of the plate (at 2 = 20 cm) as the input pressure (p in ) 

TL = 101o glo (— V (3.1) 

\PoutJ 

The present FE/BE method results are shown in Figure 3.6 and exact solutions^ 17 ! and 
experimental data^ 9 ' in Figure Note the FE/BE method results are of comparable 

accuracy even for the first mesh case. 

The second example considered analyzes the same plate/cavity system as the previous 
one. The cavity is still cubic with all five walls being acoustically rigid; however, this 
time the panel is made of composite materials. The cavity and composite plate used in 
this example have the following properties: 


Panel/Cavity Problem — Simply Supported Composite Plate 

Cavity 


X axis length 
Y axis length 
Z axis length 
Density of Air p 0 = 


L x = 8 in. 

L y — 8 in. 

L z = 8 in. 

0.1138 x 10" 6 lb/ in. 3 
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Sound speed c = 1.3504 x 10 4 in./ sec. 

Composite Plate 
X axis length = 8 in. 

Y axis length = 8 in. 

Thickness = 0.036 in. 

Graphite/Epoxy 

E\i = 2.25 x 10 7 psi 
E 22 = 1-17 x 10 6 psi 
v>i 2 = 0.22 

Gi 2 = 0.66 x 10 6 psi 
p = 0.1458 x 10~ 3 lb/ in} 

Three 6-layer composite plates of different stacking sequences are investigated. The 
plate is simply supported. The transmission losses are plotted in Figures 3.8 to 3.10 versus 
frequency. It can be seen that, in the three different stacking sequences, the fundamental 
natural frequencies increase from uncoupled 115.3 Hz for (0/30/0) s plate, 114.9 Hz for 
(0/60/0) s plate, 109.5 for (0/90/0) s plate (obtained by 7x7 mesh), to coupled 143.9 Hz, 
143.9 Hz, and 139.7 Hz, respectively. 

Boundary conditions are another factor affecting coupled natural frequencies. Figure 
3.1 1 shows the transmission losses of the system with the same (0/90/0 ) 3 laminated plate, 
but with different supported conditions. The 7x7 mesh is used. The coupled fundamental 
natural frequency increases to around 240 Hz for the clamped case. Generally, the 
transmission loss for the clamped case is larger than that of the simply supported case. 
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Figure 3.6 Transmission loss for the system consists of a simply supported isotropic 
plate and a cubic cavity: Present results 
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Figure 3.8 Transmission loss for the system consists of a simply supported (0/90/0) 
plate and a cubic cavity 




Frequency (Hz) 

Figure 3.9 Transmission loss for the system consists of a simply supported (0/60/0) 
plate and a cubic cavity 




Figure 3.10 Transmission loss for the system consists of a simply supported (0/30/0) s 
plate and a cubic cavity 
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Frequency (Hz) 


Figure 3.1 1 Transmission loss for the system consists of a simply supported (0/90/0 ) s 
plate and a cubic cavity, Comparison of different plate boundary conditions. 
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3.3 Structural-Acoustic Interaction Response by Modal Formulation 


The next objective is to form the model in modal coordinates for the coupled system 
consisting of a 6-layer (0/90/0) a composite plate and a cubic acoustic cavity, and to 
investigate the response of this coupled system under transient acoustic pressure. It is 
essential for a mathematical model to accurately represent transient behavior if the model 
is to be used as the basic equation for control. Even under harmonic inputs, there is a 
transient period immediately following controller activation. Although the final response 
of the coupled system will be harmonic under harmonic disturbance, the capability of the 
coupled equations to model the response under transient load is necessary. Meanwhile, 
to predict transient response itself is rather important for structural-acoustic interactions. 

Although the coupled finite element and boundary element formulation developed 
in Chapter 2 can be applied to transient problems directly, the current formulation is 
not convenient because of the large number of equations, and it is very difficult, if not 
impossible, to use this formulation in any control attempt. Hence, the coupled equations 
have to be transformed into the modal formulation. 

3.3.1 Coupled Natural Frequency and Mode Shape 

To transform the coupled finite element and boundary element equations in terms 
of physical coordinates into modal coordinates, a transformation matrix is required. 
As mentioned in Chapter 2, we use the eigenvectors of the coupled system as the 
transformation matrix. Generally, not all boundary element methods applied in acoustics 
can be used to obtain the natural frequency and/or modal shape of the system, because 
the frequency of the external load is embedded in the fundamental solution 1481 . But the 
dual reciprocity boundary element method is one of those which can isolate eigenvalues 
and eigenvectors. In fact, assuming the load term and damping term in equation (2.67) 
are zero, we obtain an eigenvalue problem. 

Table 3.1 gives a comparison of the natural frequency results for the coupled system 
consisting of a cubic acoustic cavity with all five walls rigid and a simply supported brass 
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plate obtained from experiment^, analytical approach 1171 , and present method. The table 
indicates the coupled FE/BE method gives rather accurate results. 


Table 3.1 Natural frequencies (Hz): Exact solution. Experiment data, and Present Results 


Exact [17] 

Experiment [9] 

Present FEM/BEM 

FE-25, BE- 125 FE-49, BE-294 

Uncoupled 

87.0 

91.0 

87.7 

87.4 

78.1 

390.4 

397.0 

394.7 

392.1 


702.5 

730.0 

710.9 



860.0 

864.0 

891.8 

872.9 

857.5 


The convergence of the time-domain modal formulation, that is how many modes 
should be retained in the analysis, was studied. Figure 3.12 shows the transmission loss 
of the system consisting of a cubic acoustic cavity with all five wall rigid and a simply 
supported brass plate under uniformly distributed external pressure. It indicates that the 
lowest five modes are sufficient to give an accurate result. For the system with a 6-layer, 
simply supported (0/90/0) a composite plate, the number of modes needed for accurate 
results is also five as shown in Figure 3.13. The non-dimensional plate deflection and 
transmitted acoustic pressure on the cavity back wall and one side wall (for example, x 
= 8 in.) for the first mode are plotted in Figures 3.14 to 3.16 to show the coupling affect 
on both the deflection mode of plate and the acoustic pressure mode of the cavity. We 
thus will take five modes for all subsequent calculations using modal formulation. 
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Transmission Loss (dB) 



Figure 3.12 Convergence of the time-domain formulation for the system consists of 
a simply supported brass plate and a cubic cavity 



Transmission Loss (dB) 
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Frequency (Hz) 


Figure 3.13 Convergence of the time-domain formulation for the system consists of 
a simply supported (0/90/0 ) 3 plate and a cubic cavity 
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a* 



Deflection 



Figure 3.14 The first modal shape of the system consisting of a simply supported 
(0/90/0) 3 plate and a cubic cavity: plate deflection 
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Acoustic Pressure 



Figure 3.15 The first modal shape of the system consisting of a simply supported 
(0/90/0) s plate and a cubic cavity: pressure on the back wall of the cavity 
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Acoustic Pressure 



Figure 3.16 The first modal shape of the system consisting of a simply supported 
(0/90/0) s plate and a cubic cavity: pressure on a side wall of the cavity 
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3.3.2 Time-Domain Modal Results 


When modal coordinates, or modes, are selected, we are ready to predict the response 
of the coupled structural-acoustic system. The time-domain modal formulation is actually 
a set of ordinary differential equations (for transient problem). Since the number of 
equations in the modal formulation is small (here it is five), a classical Runge-Kutta 
numerical integration scheme can be easily employed. 

What we study here is the coupled structural-acoustic system consisting of a 6-layer, 
(0/90/0) 3 composite plate and a cubic cavity given in Section 3.2 under uniformly 
distributed pressure on the external surface of the plate. The initial conditions in all three 
examples are static. The first example problem is the response to an impulse, that is, 
the forcing function is assumed to be 

F(x,y,t) = 1.56(f) x 10 -3 ( Ib/in . 2 ) (3.2) 

The time history of the displacement at the center of the plate, the transmitted acoustic 
pressure at the center of the cavity back wall (x = y = 4 in. z = 0), and the forcing 
function are shown in Figure 3.17. It can be seen that multiple modes are excited. 

Multiple resonant excitation of this coupled structural — acoustic system is studied 
next. The forcing function is 


f{x,V,t) = 1.5[sin (3007 t<) + sin (1760 tt<)] x 10~ 3 ( Ib/in , 2 ) (3.3) 

and, hence, it excites the first and fourth system modes which have natural frequencies of 
139.7 Hz and 879.0 Hz, respectively. This forcing function represents a periodic plane 
wave with a sound pressure level of 126 dB. The time histories of the forcing function, 
the displacement at the center of the plate, and the transmitted acoustic pressure at the 
center of the cavity back wall are shown in Figure 3.18. 


The last example of this chapter is the response under a forcing function defined as 
f 95 1 0 < f < 0.01 


f(x,y,t) = 1.5 x ( 


0.95 

0.95 - 55(f - 0.03) 
0.4 


0.01 < t < 0.03 
0.03 < f < 0.04 

t > 0.04 


> (Ib/in. 2 ) (3.4) 
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The forcing function, the displacement at the center of the plate and the transmitted 
acoustic pressure at the center of the cavity back wall are shown in Figure 3.19. Those 
examples demonstrated the capability of handling time-dependent problems of the present 
method. 



Time (sec.) 


Figure 3.17 Displacement at the center of the plate and transmitted pressure at the 
center of the cavity back wall: Forcing function (bottom) equation (3.2) 
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Pressure (lb/in. A 2) Displacement (in.) 



Figure 3.18 Displacement at the center of the plate and transmitted pressure at the 
center of the cavity back wall: Forcing function (bottom) equation (3.3) 
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Pressure (lb/in. A 2) Displacement (in.) 



Time (sec.) 


Figure 3.19 Displacement at the center of the plate and transmitted pressure at the 
center of the cavity back wall: Forcing function (bottom) equation (3.4) 
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3.4 Conclusions 


The numerical examples given in this chapter show that the dual reciprocity boundary 
element method, the coupled boundary element-finite element method and its time-domain 
formulation, developed in Chapter 2, can be used to solve acoustic or coupled structural- 
acoustic problems accurately and efficiently. Although we did not give any example 
with complex geometry, based on the versatility of the boundary element method and 
the finite element method, and the general approach employed in this chapter for the 
structural-acoustic coupling, we can expect that it is not difficult to extend this method 
to problems of complex geometry. In fact, the example problems investigated in this 
chapter show little difficulty in applying this method to different plates, different plate 
boundary conditions, and different acoustic boundary conditions. 
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Chapter 4 

CONTROL OF INTERIOR ACOUSTIC PRESSURE LEVEL 


The second step of this investigation is to control the interior acoustic pressure level 
based on the governing equations in the time-domain derived in Chapter 2. The control 
force is applied through the piezoelectric patches bonded on the two surfaces of the plate. 

The governing equations in time-domain are first rewritten in state space form. Then 
a feedback controller, the linear quadratic regulator, and a feedforward controller based on 
the concept of least mean square (LMS) are employed to reduce the acoustic pressure field 
inside the cavity via transmission through the flexible plate using piezoelectric patches 
as actuators. 


4.1 State Space Formulation 

The time-domain modal formulation is rewritten here as 

[J/]*M + KTW + [ATM = {ii} + {M (4.D 

where { L \ } is due to external acoustic pressure and {£ 2 } is due to piezoelectric actuators. 
If we denote the system response due to external pressure as {/?«} and that due to actuators 
as {r) a }, then we have 

WM} + m-i,} + ffw = (4.2) 

[M]*{rja) + (CH"?.} + [ATM) = {Li} (4.3) 

and 

{*)} = M + {>?<»} ( 4 - 4 > 

Define 

Xo = [»?e,7?e] r , X = [7/a,f/a] T ( 4 - 5 ) 

Then we have the state space form of the governing equations (4.3) as 

x = Ax + Bu (4.6) 
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where u — {E%} / E 2max is the control variable vector, Ezmax is the maximum allowable 
operating electric field of the piezoelectric material, 

A= -M'~ X C* ,B= (47) 

and [G] is the control inference matrix which is obtained from 

{L 2 } = [G\ u (4.8) 

4.2 Linear Optimal Controller 

The linear quadratic regulator is considered in this section to reduce the acoustic 
pressure level inside the cavity. The objective here is to minimize the global acoustic 
pressure using the piezoelectric patches bonded on the two surfaces of the plate. The 
control mathematical models based on the time-domain state space formulation are 
presented first, then numerical results are obtained. 

4.2.1 Mathematical model 

Define the linear quadratic performance index for optimal control as^O 

OO 

J=\j (x + x 0 ) T Q(x + Xfl) + u t Ru (ft (4.9) 

o 

where Q is a real symmetric positive semi-definite matrix and R is a real symmetric 
positive definite matrix. Then the optimal control for this linear quadratic problem can 
be obtained as a function of the costate^ 8 ^, and if we assume the relation between the 
costate and the state variable is linear, we have the control effort as 

u = -R -1 B t P(x + x 0 ) (4.10) 

while the Riccati matrix P, which is symmetric and positive definite, can be obtained 
from the Riccati equation^ 82 ! 

A t P + PA - PBR _1 B t P = -Q (4.11) 
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Consider a maximum electric field that can be applied to the piezoelectric material, 
where the constraint for equation (4.11) is 


Ui| < 1, i — 1*2, Np 


(4.12) 


and Np is the number of piezoelectric actuators. Then equations (4.9), (4.10), and (4. 1 1) 
form a bounded optimization problem. 

Since our objective is to control the inside acoustic pressure level, we define 
the weighting matrix Q so that the (x + Xo) T Q(x + x c ) term in the linear quadratic 
performance index represents the acoustic potential energy in the cavity 

( 4 . 13 ) 

V 

where ['Kn] is the modal shape matrix evaluated inside the cavity, it relates to the NN 
interior points in the cavity and can be obtained from the coupled modal matrix ['P ] in 
equation (2.70). Then we have 


[Q] = 


Qn 

o 


o 

Q22 


(4.14) 


where 1 

Qu = 

2p 0 c i 

Q22 = 0 


(4.15) 


4.2.2 Numerical Results 


Based on the above mathematical model, several example problems are solved. The 
coupled structural-acoustic system is the cubic cavity backed by a simply supported 
6-layer (0/90/0 ) 3 composite plate given in Section 3.2. The results are presented as 
follows. 


Piezoelectric Layers 

Before the discussion of control results, the coupled structural-acoustic system is 
modified by bonding a piezoelectric layer on each of the external and internal surfaces 
completely of the composite plate. The piezoelectric material is considered as isotropic. 
The characteristics of the piezoceramic layer are: 1711 
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Piezoelectric Layer (PZT) 

Modulus of Elasticity E p = 0.9 x 10 7 psi 
Poisson's Ratio j p = 0.3 
Density p p = 0.7101 x 10 -3 lb/in. 3 

Stress / Charge Constant d n = d 32 = -7.51 x 10 -9 in./v 
Max. Electric Field E 3max = 1.52 x 10 4 v/in. 

The thickness of the piezoelectric layer is assumed to be the same as the composite 
layer; that is, 0.006 inch. Following the same procedure given in Chapter 3, the 
responses of the full coupled acoustic/composite/piezoelectric system under harmonic 
acoustic pressure with different frequency content are obtained. The transmission loss at 
the center of the cavity back wall is compared with that of the system without piezoelectric 
layers in Figure 4. 1 . The natural frequencies of the systems with or without piezoelectric 
layers are tabulated in Table 4.1. Two conclusions can be obtained from this result: 1) the 
mass and stiffness of the piezoceramic layers have to be considered, the responses of the 
structural-acoustic systems with and without the piezoceramic layers are rather different, 
the natural frequencies of the system are changed, the fundamental frequency decreases 
10Hz and the second increases because of the isotropic property of the piezoelectric 
layers; and 2) five modes are sufficient for the time-domain modal formulation to obtain 
accurate results for input frequency up to 1000 Hz. 


Table 4.1 Natural frequencies (Hz) of the system with and without piezoelectric layers 



mode 1 

mode 2 

mode 3 

mode 4 

mode 5 

With PZT Layers 

129.6 

574.7 

685.0 

867.9 

1060.1 

Without PZT Layers 

139.7 

509.8 

783.5 

878.6 

985.4 


Resonant Excitation 

The forcing function for this example is assumed to be 

f(x,y,t ) = 3.0 sin (2807ft) x 10 -3 (lb/in 2 .) (4.16) 
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which corresponds to the frequency of 140 Hz and, hence, the fundamental natural 
frequency of the coupled structural-acoustic system with piezocermatic layers which is 
129.6 Hz is excited. Consider the reflection, equation (4.16) models a periodic exterior 
plane wave with a maximum pressure level of 120 dB. The initial condition of the 
system is assumed to be at rest. The weighting matrix R is taken to be diagonal, and 
each diagonal element is equal to 10 -4 . 

Table 4.2 Modal coefficients for controlled and uncontrolled system for resonant 
excitation 



mode 1 

mode 2 

mode 3 

mode 4 

mode 5 

uncontrolled 

1.77E-01 

1.27E-04 

1.87E-04 

-6.34E-04 

-2.34E-05 

controlled 

3.91E-03 

4.56E-04 

5.29E-04 

-5.11E-04 

-8.56E-05 


Control is implemented via piezoelectric patches covering the two surfaces of the 
composite plate. The time histories of the displacement at the center of the plate, the 
pressure at the center of the cavity back wall, the norm of the pressure field inside the 
cavity (||[ 1 fn]{f?}||), and the control effort are given in Figure 4.2. The amplitude of the 
acoustic pressure at the center of the cavity back wall decreases from 7.22 x 10 -3 lb jin. 2 
to 1.37 x 10~ 4 lb jin. 2 after the controller is turned on, that means a 34.4 dB increase in 
the transmission loss. At the same time, the amplitude of the maximum deflection of the 
plate reduced from 7.93 x 10 -4 in. to 1.14 x 10 -5 in. It can be seen form Figure 4.2, 
the acoustic pressure and the displacement of the plate are all reduced when the control 
force is applied. Figures 4.3 and 4.4 give the comparison of the maximum acoustic 
pressure inside the cavity and the maximum plate displacement between the controlled 
and uncontrolled system at the instant of maximum acoustic pressure at the center of 
the cavity back wall. It can be seen that, for the uncontrolled system, the response 
is dominated by the first coupled mode. After the controller is turned on, multi- mode 
behavior appears. Table 4.2 gives the values of the modal variables corresponding to 
maximum acoustic pressure at the center of the cavity back wall for the controlled and 
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uncontrolled systems, those results indicate the dominated fundamental mode is reduced 
with active control. 

Multi- Resonant Excitation 

The forcing function for this example is 

f(x,y,t) = 3.0[sin (280 tt/) + sin (1760x0] x 10" 3 (/6/m. 2 ) (4.17) 

Hence, the fundamental and the fourth coupled natural frequency, the later corresponds 
to the uncoupled cavity frequency at 867.9 Hz, are excited. This function models a 
periodic plane wave of maximum pressure level around 126 dB. The initial condition of 
the coupled system is at rest. 

Table 4.3 Modal coefficients for controlled and uncontrolled system for multi-resonant 
excitation 



mode 1 

mode 2 

mode 3 

mode 4 

mode 5 

uncontrolled 

1.78E-01 

3.68E-04 

5.48E-04 

-4.43E-02 

-6.48E-05 

controlled 

4.01E-03 

6.15E-04 

9.10E-04 

-2.82E-03 

-4.61E-05 


The same cost function as in the resonant excitation case is applied in this example. 
The time histories of the displacement at the center of the plate, the pressure at the center 
of the cavity back wall, the norm of the pressure field inside the cavity, and the control 
effort are given in Figure 4.5. Notice that the pressure and displacement are both reduced 
dramatically after the controller is turned on. The amplitude of the acoustic pressure at the 
center of the cavity back wall decreases from 9.56 x 10 -3 lb/ in. 2 to 2.33 x 10 -4 lb/ in. 2 
after the controller was turned on, that means a 32.3 dB increase in the transmission loss, 
and the amplitude of the maximum deflection of the plate reduced from 8.05 x 10 -4 in. 
to 1.50 x 10 -5 in. Figures 4.6 and 4.7 give the comparison of acoustic pressure inside 
the cavity and plate displacement of the controlled and uncontrolled system at the instant 
of maximum acoustic pressure at the center of the cavity back wall. Table 4.3 tabulates 
the modal coefficients corresponding to maximum acoustic pressure at the center of the 
cavity back wall for the controlled and uncontrolled systems. It can be seen that the 
fundamental and fourth modes are the dominate and both are reduced. 
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Off-Resonant Excitation 


In this example, the forcing function was taken to be 

/(ar, y, t) = 3.0 sin (5007r< ) x 10 -3 (/6/m. 2 ) (4.18) 

with a frequency of 250 Hz which is between the natural frequencies of the first (129.6 
Hz) and the second (574.7 Hz) system modes. As noted before, this models a periodic 
plane wave of maximum pressure level around 120 dB. The initial condition of the 
system is again at rest. 

The same quadratic cost as before is applied in this example. The time histories of the 
displacement at the center of the plate, the pressure at the center of the cavity back wall, 
the norm of the pressure field inside the cavity, and the control effort are given in Figure 
4.8. Figures 4.9 and 4.10 give the comparison of acoustic pressure inside the cavity and 
the plate displacement of the controlled and uncontrolled system. The amplitude of the 
acoustic pressure at the center of the cavity bottom decreases from 6.49 x 10 -4 /6/m. 2 
to 1.40 x 10 -4 /6/m. 2 after the controller was turned on, that means a 13.3 dB increase in 
the transmission loss, and the amplitude of the maximum deflection of the plate reduced 
from 6.49 x 10~ 5 in. to 1.92 x 10 -5 in. Notice that the reduction of acoustic pressure 
and the plate deflection in this off-resonant excitation example is much smaller than 
in the resonant or multi-resonant examples, this phenomena is reasonable because, the 
pressure or deflection of the off-resonant case themselves are much small comparing to 
resonant or multi-resonant cases. This can be observed by comparing Table 4.4 with 
Tables 4.2 and 4.3. 

Table 4.4 Modal coefficients for controlled and uncontrolled system for off-resonant 
excitation 



mode 1 

mode 2 

mode 3 

mode 4 

mode 5 

uncontrolled 

1.52E-02 

1.54E-04 

2.02E-04 

-9.13E-04 

-2.36E-05 

controlled 

3.46E-03 

4.01E-04 

4.40E-04 

-7.86E-04 

-5.89E-05 
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Effect of Weighting Matrix Q 


Weighting matrix Q determines the linear quadratic response term in the cost function. 
Consequently, different definitions of Q will give different performance. Here we defined 
another weighting matrix to compare the performance with the cavity acoustic potential 
energy given in equation (4.15). The elements in equation (4.14) are defined here as 


Qu = [K]* 
Q 22 = [M]* 


(4.19) 


The results of the system response under external excitations given in equations 
(4.16) to (4.18) are obtained using the same R matrix as before. The time histories of 
the displacement at the center of the plate, the pressure at the center of the cavity back 
wall, the norm of the pressure field inside the cavity, and the control effort are given in 
Figures 4.11 to 4. 1 3 for the three loads, respectively. Table 4.5 compares the increase 
of sound transmission loss at the center of the cavity back wall and the displacement 
reduction at the center of the plate for those two Q matrices. The displacement reduction 
is defined as the ratio of the controlled and the uncontrolled displacement amplitudes. It 
can be seen that the performance from the two different cost functions are comparable 
to each other. 


Table 4.5 The increase of sound transmission loss (ITL) at the center of the cavity 
back wall and the displacement reduction (DR) at the center of the plate: comparison 
of different weighting matrices Q 



Q matrix in Eq. (4.15) 

Q matrix in 

Eq. (4.19) 

ITL (dB) 

DR (%) 

ITL (dB) 

DR (%) 

Resonant 

34.4 

1.44 

27.8 

4.04 

Multi-Resonant 

32.3 

1.86 

29.2 

4.35 

Off-Resonant 

13.3 

29.4 

12.4 

29.3 
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Figure 4.1 Comparison of responses to harmonic external excitations for the 
coupled structural-acoustic system with or without piezoelectric layers. 
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Figure 4.2 Time history of the pressure at the center of the cavity back wall, 
the displacement at the center of the plate, the norm of the cavity pressure field, 
and the control effort: Resonant excitation 
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Figure 4.3 Uncontrolled pressure field and plate deflection: Resonant excitation 
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Figure 4.4 Controlled pressure field and plate deflection: Resonant excitation 
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Figure 4.5 Time history of the pressure at the center of the cavity back wall, 
the displacement at the center of the plate, the norm of the cavity pressure field, 
and the control effort: Multi-Resonant excitation 
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Figure 4.6 Uncontrolled pressure field and plate deflection: Multi-Resonant excitation 
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Figure 4.10 Controlled pressure field and plate deflection: Off-Resonant excitation 
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Figure 4.1 1 Time history of the pressure at the center of the cavity back wall, 
the displacement at the center of the plate, the norm of the cavity pressure field, 
and the control effort: Resonant excitation and Q matrix in equation (4.19) 
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Figure 4.12 Time history of the pressure at the center of the cavity back wall, 
the displacement at the center of the plate, the norm of the cavity pressure field, 
and the control effort: Multi-Resonant excitation and Q matrix in equation (4.19) 
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Figure 4.13 Time history of the pressure at the center of the cavity back wall, 
the displacement at the center of the plate, the norm of the cavity pressure field, 
and the control effort: Off-Resonant excitation and Q matrix in equation (4.19) 
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4.2 3 Optimal Location of the Piezoelectric Patch 

We are now trying to investigate the best location of the piezoelectric patch. In 
this study, we will utilize the finite element method to determine the locations where 
piezoelectric patches can most affect the inside acoustic pressure. To assess the effect 
of the different patch locations, a norm of the feedback control gain (NFCG) is defined. 
The feedback control gain matrix can be determined from equation (4.10) as* 71 l 

[G.] = [R\~'\B] T [P\ (4.20) 

the norm of the feedback control gain is defined from the norm of matrix [G a \ as 

Np 2m 

kfcg=J: 

<=i L> =i 

where Np is the number of piezoelectric actuators, and 2m is the total number of state 
variables in equation (4.6). 

Assuming that we use a pair of piezoelectric patches (Np = 2) having the same shape 
and area of a finite element and bonded to the two surfaces of the composite plate. The 
NFCGs for the piezoelectric patches covering each element alone are calculated. The 
same weighting matrices Q and R as previously defined are used. After normalization 
with the maximum norm to be unity, they are plotted in Figure 4.14. Only 16 elements 
are marked and the others can be obtained by symmetry. It can be seen that placing the 
actuators at the center of the plate will give the best performance. 

The results of using a pair of actuators at the center of the plate for the three external 
loads are given in Figures 4.15 to 4.17, respectively, and, as shown in Table 4.6, they are 
compared with the previous results in the sense of using the same weighting matrices. 
Notice that the optimally located actuators uses only 2% of the piezoelectric material 
as the earlier analysis and computations with piezoelectric layers covered the complete 
surfaces. 


(4.21) 
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Table 4.6 The increase of sound transmission loss (ITL) at the center of the cavity 
back wall and the displacement reduction (DR) at the center of the plate: 
Comparison of the performance of optimal located actuators and previous ones. 



Actuators cover two 
surfaces of plate 

Actuators cover two 
surfaces of center element 

ITL (dB) 

DR (%) 

ITL (dB) 

DR (%) 

Resonant 

34.4 

1.44 

18.1 

13.3 

Multi-Resonant 

32.3 

1.86 

23.1 

13.9 


13.3 

29.4 

8.70 

40.0 


0.09 

0.18 

0.06 

0.31 




0.23 

0.26 

0.13 

0.40 




0.23 

0.23 

0.21 

0.61 




0.53 

0.55 

0.41 

1.00 


























Figure 4.14 NFCGs for different locations of piezoelectric actuator. 
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Figure 4.16 Time history of the pressure at the center of the cavity back wall, 
the displacement at the center of the plate, the norm of the cavity pressure field, 
and the control effort: Multi-Resonant excitation and best location 
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Figure 4.17 Time history of the pressure at the center of the cavity back wall, 
the displacement at the center of the plate, the norm of the cavity pressure field, 
and the control effort: Off-Resonant excitation and best location 
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4.3 Feedforward Controller 


In linear quadratic regulator, the control gain is sought to give the best trade-off 
between performance level and cost of control. However, this controller requires full 
state feedback which in reality is hard to achieve^ 83 !. In this section, a more practical 
controller design, feedforward control^ 58 !, is investigated. 

Assume the external acoustic wave is measurable, we try neutralize the external 
disturbance by applying a control force (from piezoelectric actuators) in the weighted 
least-mean-square sense which counteracts the external pressure disturbance. 

4.3.1 Mathematical Model 


Assume that the external wave is a plane wave, and the measured pressure on the 
plate can be expressed as a Fourier series 

Nl 

P(x, y,t) = ^2 P k sinwjfct (4.22) 

k=\ 

where Pk are known coefficients, and u>k are the input frequencies. Then the load term 
in equation (4.2) can be expressed as 

N l 

(M = {M 52 Pk sinwjfct (4.23) 

i t=i 

where {Zi } is the modal force corresponding to uniformly distributed unit pressure. 

To suppress the acoustic field inside the cavity, we assume that the control force 

has the form of 

Nl 

{L 2 } = [G]{/3} 52 sinw ifc< ( 4 -24) 

k=i 


where [G] is the modal force matrix corresponding to the maximum operating electrical 
field as defined in equation (4.8). Each column of [G] corresponds to each actuator. 
Further, {/?} is the undetermined coefficient vector. 

After substituting equation (4.23) and (4.24) into equation (4.1), we have 

N l 

\M]’W + m-i} + [ATM = ({£ 1 } + [G]{/?}) £Pisinu, t l (4.25) 

k=\ 
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If we ignore the damping term when we determine the control effort, and notice that 
matrices [ M ]* and [ K ]* are all diagonal, we have 

N l N l 

M = £ {a} fc sinu^i = + [G]{0}) sinu k t (4.26) 

h=l k=l 

where matrix [0] fc is also diagonal and 


Mu = 


Pk 


k*ii ~ ^l m U 


(4.27) 


To reduce the acoustic pressure level, we desire the inside acoustic pressure caused 
by external acoustic waves at each frequency to be null, that is 


['J'nKahsinwjt* = {0}, k = 1,2,..., N L (4.28) 


Generally, the number of the actuators is not the same as the number of frequencies in 
the exciting wave, and there is no proportional relation between {L\} and [G]. Thus, 
equation (4.28) is impossible to achieve. Instead, {/?} will be designed to reduce the 
internal pressure waves as much as possible in the least mean square sense, or 

N l N l 

n = £ ll[i'niM i ll 2 = £ ||[®n][0U({M + [GH«)|| 2 - mm. (4.29) 

*=1 Jfc=l 

The resulting {0} is 

E N T [0][[* n ] T [» n i[0] 1 

m = ~wr f Zi > (4 - 30) 

E [Gl T [0]r[f n ] T ['ffi][0] 1 [GI 

k=l 

4.3.2 Numerical Results 

The example problems in Section 4.2 using three different excitation forces investi- 
gated are studied again here using the feedforward controller. As in the LQR, only one 
pair of actuators covering the external and internal surfaces of the composite plate are 
used. The results are presented as follows. 
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Resonant Excitation 


The comparison of the controlled and uncontrolled responses is given in Figures 
4.18 and 4.19. The feedforward controller reduced the maximum pressure at the center 
of the cavity back wall from the uncontrolled 7.22 x 10 -3 /6/m. 2 to controlled 

5.69 x 10 -5 /6/m. 2 This is a 42.1 dB increase in the transmission loss. At the same time, 
the amplitude of the maximum deflection of the plate is reduced from 7.93 x 10 -4 in. 
to 2.05 x 10 -5 in. Higher performance has been achieved here relative to the LQR case. 

Multi-Resonant Excitation 

The results are compared in Figures 4.20 and 4.21. The feedforward controller 
reduced the maximum pressure at the center of the cavity back wall from the uncontrolled 
9.56 x 10 -3 /6/m. 2 to controlled 2.62 x 10 -3 /6/m. 2 This is a 11.2 dB increase in the 
transmission loss. The amplitude of the maximum deflection of the plate is reduced from 
8.05 x 10" 4 in. to 1.23 x 10 -4 in. 

Off-Resonant 

The results are given in Figures 4.22 and 4.23. The feedforward controller reduced 
the maximum pressure at the center of the cavity back wall from the uncontrolled 
6.49 x 10 -4 lb/in. 2 to controlled 1.17 x 10 -4 /6/m. 2 This is a 14.9 dB increase in 
the transmission loss. At the same time, the amplitude of the plate deflection is reduced 
from 6.51 x 10 -5 in. to 3.34 x 10 -5 in. 
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Figure 4.18 Comparison of controlled (left) and uncontrolled (right) system response 
to resonant excitation: Time history of the pressure at the center of the cavity back 
and displacement at the center of the plate 
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Figure 4.20 Comparison of controlled (left) and uncontrolled (right) system response 
to multi-resonant excitation: Time history of the pressure at the center of the cavity 
back wall and displacement at the center of the plate 
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Figure 4.21 Comparison of controlled (left) and uncontrolled (right) system response 
to multi-resonant excitation: norm of the cavity pressure field, and control effort 
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Figure 4.22 Comparison of controlled (left) and uncontrolled (right) system response 
to off-resonant excitation: Time history of the pressure at the center of the cavity 
back wall and the displacement at the center of the plate 
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Figure 4.23 Comparison of controlled (left) and uncontrolled (right) system response 
to off-resonant excitation: the norm of the cavity pressure field and the control effort 
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433 Optimal Location of the Piezoelectric Patch 

Once again, we try to find the optimal location of the piezoelectric patches. As we 
mentioned in Section 4.2, the mass and stiffness of the piezoceramic patches will affect 
the characteristics of the coupled system. Thus, the different location and shape of the 
piezoelectric actuators will mathematically change the system and the response. So in 
this study, we only try to investigate the optimal location when the external excitation is 
a simple harmonic plane wave at frequency of 140 Hz. And similar to Section 4.2, we 
try only one pair of piezoelectric patches bonded to the two surfaces of the plate, and 
each patch occupies only one finite element. 

The criteria of the performance in this section is defined as the increase of the 
transmission loss for the norm of the inside acoustic pressure field between controlled 
and uncontrolled systems. The results are given in Figure 4.24 after normalization with 
the maximum value. Notice the best location of the piezoceramtic actuator is at the center 
of the plate again. The time histories of the displacement at the center of the plate, the 
pressure at the center of the cavity back wall, the norm of the pressure field inside the 
cavity, and the control effort when the actuator is at the best location are given in Figures 
4.25. The comparison of the increase of sound transmission loss at the center of the 
cavity back wall and displacement reduction at the center of the plate between actuators 
at the best location and actuators covering two surfaces of the plate are given in Table 4.7. 

Table 4.7 The increase of sound transmission loss (ITL) at the center of the cavity back 
wall and the displacement reduction (DR) at the center of the plate: Comparison of the 
performance of optimal located actuators and the one covering two surface of the plate. 


Actuators cover two surfaces of plate 

Actuators cover center element 

ITL (dB) 

DR (%) 

ITL (dB) 

DR (%) 

42.1 

2.59 

33.8 

5.37 
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Figure 4.24 The normalized increase of transmission loss obtained from different location 
of piezoelectric actuators. 
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Figure 4.25 Comparison of controlled (left) and uncontrolled (right) system responses 
to resonant excitation: Time history of the pressure at the center of the cavity back 
and displacement at the center of the plate: Best location 
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4.4 Conclusions 


The linear quadratic regulator and a feedforward controller based on the concept of 
least mean square are applied with the intent to reduce the acoustic pressure field inside 
the cavity for the coupled plate/cavity system under different external disturbances. The 
piezoelectric patches bonded on the plate surfaces are used as the actuators, and the 
total acoustic potential energy inside the cavity is used as the control objective. The 
pressure level inside the cavity and the plate deflection are reduced successfully for both 
controllers. The optimal location of the piezoelectric actuators are investigated based on 
those two controllers. 
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Chapter 5 

Summary and Conclusions 
5.1 Concluding Remarks 

The contributions of this dissertation consists of two parts, one is the modeling of 
the acoustic-structure interaction using the coupled finite element and boundary element 
methods, the other is the suppression of interior acoustic pressure levels transmitted 
through composite plates using bonded piezoelectric actuators. 

To predict the response of the coupled structural-acoustic system under external 
acoustic excitation, the C 1 conforming finite element method is used to model the plate, 
and the dual reciprocity boundary element method is used to model the interior acoustic 
domain. Based on the continuity requirement on the interior surface of the plate, the 
boundary element method and finite element method are coupled together to form a 
powerful tool which can be used to solve structural-acoustic interaction problems. Interior 
acoustic problems are analyzed first using only the boundary element method, serving to 
verify that the boundary element method is an accurate and versatile approach in solving 
acoustic problems. The coupled method is then applied to a cubic acoustic cavity backed 
by a brass plate, the results are compared with the analytical solution and experiment 
data, and the accuracy of the coupled method is verified. This method is then applied to 
the cavity backed by a composite plate. 

Taking advantage of the dual reciprocity boundary element method, the coupled finite 
element and boundary element scheme is used to calculate the eigenvalues and eigen- 
vectors of the coupled structural-acoustic system. Those eigenvalues are the natural 
frequencies of the coupled system. The eigenvectors are employed to form the transfor- 
mation matrix to transform the coupled equations into modal formulation which has a 
very small number of equations compared to the original one. The responses of the cou- 
pled acoustic-structure system under various external excitations are investigated using 
the modal formulation to demonstrate its accuracy and capability. 
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The results show that it is necessary to consider the interaction between the acoustic 
domain and the structure, and also to consider the mass and stiffness of the piezoelectric 
patches. The interaction between structure and acoustic cavity significantly influences 
the overall system natural frequencies, and the full coupling between composite laminate 
and piezoelectric patches also modifies the natural frequencies of the system. 

Based on the modal formulation, two different control designs are applied to reduce 
the transmitted acoustic pressure inside the cavity with embedded piezoelectric patches 
as actuators. The linear quadratic regulator gives the best control performance under 
certain information of the system and external load. The feedforward controller is based 
on the attempt to neutralize the external excitation in the weighted least mean square 
sense. For the two controllers, they both can reduce the acoustic pressure field inside 
the cavity and the deflection of the plate effectively. The higher the transmitted pressure 
is, the better the control performs. Specifically, for the resonant excitation, the LQR 
can increases the transmission loss up to 34.4 dB, and feedforward, up to 42.1 dB, at 
the same time, the maximum deflection of the plate is reduced 98.6% using LQR, and 
97.4% for feedforward controller. 

The optimal locations of the piezoelectric actuators are investigated for those two 
controllers. Comparison to the case where piezoelectric patches cover the entire external 
and internal surfaces of the plate, the optimal located piezoelectric actuators perform 
rather efficiently. 

Currently, industry is striving to reduce the interior noise in aerospace and automobile 
systems. These goals require more efficient approaches to model the interaction between 
acoustic and structures, and also the design of suitable controllers. This study is only an 
attempt to enhance our understanding in that direction. 

5.2 Future Work 

Many areas of research still need to be considered in order to more accurately 
model the real world systems. Present finite element/boundary element model could 
be improved in at least three areas. First, the nonlinear large displacement-strain relation 
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of the composite laminate needs to be included when external acoustic pressure is high 
or the exciting frequency is around the natural frequencies of the system. For aerospace 
engineering, the acoustic pressure level for new supersonic aircraft will be of 190 dB, the 
thin composite plate will then be into large deflection range. Second, the thermal effect 
on composite material and piezoelectric material should be considered. For supersonic 
aircraft, the surface temperature of the fuselage will easily reach few hundreds degree. 
Third, other elements, such as linear or other higher order boundary elements, or three 
dimensional elements, like shell elements, should be included to model real system as 
aircraft fuselage or automobile bodies. 

For control, other controllers should be considered. In this study, piezoelectric 
patches are used only as actuators, the use of piezoelectric patches as sensor and actuator 
becomes naturally the next step. 
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APPENDIX A 

Derivation of Finite Element Equation 


The finite element equation of motion for the laminated composite plate with piezoce- 
ramic layers can be obtained from the generalized Hamilton’ s principle. When the plate 
is discretized by certain type of elements, this principle is valid for each element. Con- 
sidering a certain element, substituting equations (2.18) to (2.21) into equation (2.17), 
we have 


j ^pSiviv -I- {6E} T {D} - {<5e} T {cr} + dV+ 
v e 

J 8wF a dS + J w8F s dS - J 6<t>qdS + 8wF c = 0 

S e i Sc3 


(A.1) 


where V e is the volume of the element, s e i is the element boundary with prescribed 
forces, s e 2 is the element boundary with given displacement, and s e3 is the surface of 
piezoelectric material. For thin composite panel, equation (A.l) becomes 

A/ 2 

pSww + {8E} T {D} — {6e} {<7 

—A/2 A < Ae Ae3 

where h is the plate thickness and p is the external loading on the plate surface A e . 



}j dAdz + J SwpdA — J 8<f>qdS — 0 (A. 2) 


The displacement vector in the element can be expressed as the displacements at 
element nodes as (see Appendix B) 



w 

■[«»]{“) 1 

' \H w )[T b }{w b } ' 

w = < 

u > = < 

[H.W 1 = ' 

[H a }[T m }{w m } 


U 1 

Jff.HM J 



and the strain— displacement relation in equation (2.22) yields the strain vector as 


{e} = [C m ]{6} + ztCfcKa} 


(A.4) 


then the constitutive equation ( 2.14) gives the stress vector for k-th layer as 
Mfc = [Q} k ([Cm]{b} + *[<?*]{«} - [d\i E *}) 
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where 


H 


k 


dxk^ik # * * d x k$jk 

dyk^ik “* d yk 6jk 

dxyk^ik * * * d X y k 8j k 


and {£ 3 } = [(£ 3 ) nl ,(£ 3 ) n2 ,- 
the location of piezoceramic. 


iT 


“1 {^3 )nn<j> 


, we assume the layers nl, n2, 


(A. 6 ) 


nn<f> are 


Substitute equations (A.4) to (A. 6 ) into equation (A.2), we have the governing 
equation in element level as ^ 711 




_ 

f » \ 

IVb 



_ 

f • 

Wb 

mf, 

0 

0 



1 

?r- 

o* 

kbm 




0 

m m 

0 

< 

Wm 

► + 

kmb 

k m 

km<f> 

< 

W m 

0 

0 

0 


0 



k(f)m 

— 1 

-se 









l"V J 


where the element matrices are defined as 



(A.7) 


K] = (p4hi + p e h e )[T b ] T J { H w }{H w fdA[T b ] (A.8) 

A e 

[m m ] = (p+ht + p c h c )[T m ff [H u ,H v ][H u ,H v ] T dA[T m ] (A.9) 

Ae 

W = Pi) T J [Cjf ([B] c + [£>],) [Cb]dA[Tb] (A. 10) 

[M = \T m f J [C m ] I '([A] e + [A\^[C m \dA[T m \ (A. 11) 

A e 

pui = ih m ] T = [r m ] T J [c m \ T ([B] e + [£»]*) nmn] <A.i2> 

Ae 

M = \B k \ T [G t \ T j [C 6 ]<M(7i] (A.13) 

Ae 

[A*»] = [B»] r j [C m )iA[T m ] (A. 14) 

Ae 

[4] = (A.15) 

{h) = {Tk] T j ?A»w\ T dA (A. 16) 
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where 


{/m} = [r m ] T J (p x [H u ) T + Py [H v ] T )d A 

At 

(A. 17) 

{U} = J <i dA 

A e 

(A. 18) 

h/2 

A 


( [-^V] » \PA ) = J [^]ki d \k^ ,Z ^ Z 

(A. 19) 


where [Bh] has been defined in equation (2.24). 

It is noticed that, in the element mass matrices (A.8 and A.9) and stiffness matrices 
(A. 10 and A.ll), the modification of piezoceramic layers is included, and p c are the 
mass density of piezoelectric material and composite material, respectively, h $ and h c 
are the thickness of piezoceramic layers and composite layers, respectively. [A], [Z?] and 
[ D] matrices for composite and piezoelectric lamina are defined as 


([A] c ,[B] c ,[£>] e ) = J [Q c ) t (l,z,z 2 )dz 
where the integrals are through the thickness of composite layers, and 

where the integrals are through the thickness of piezoelectric layers. 


(A.20) 


(A. 21) 
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APPENDIX B 

Finite Element Characteristics 


The rectangular element under consideration is with 24 structural degrees of freedom 
(DOF) for bending and membrane plus one extra electrical DOF for each piezoceramic 
layer. The unknown displacements, w , w x , w y , w xy , u, and v are defined at each 
element node, while the electric potential <t> (electrical DOF) is defined for each piezo- 
electric layer (see Figure 2.2). 


B.l Transformation between nodal displacements and generalized coordinates 


The transformation matrices between nodal displacements and generalized coordi- 
nates can be derived by considering an element in its local coordinates. For the bending 
deflection, the C l conforming element assumes the bending deflection in the element is 
distributed as: 

w(x,y, t) = a\ + a 2 x + a 3 y + a 4 x 2 + a 5 xy + a 6 y 2 + a 7 x 3 
+ a 8 x 2 y + a 9 xy 2 + a i 0 y 3 + a n x 2 y + a u xy 3 

(B.l) 

+ a\3X 2 y 2 + a u x 2 y 2 + ai 5 x 2 y 3 + a^x z y z 


= [ H w {x,y)]{a(t )} 

where [H w (x,y)} is the transverse shape function 

[H w {x,y)\ = [l,x,y,x 2 ,xy,y 2 ,x 3 ,x 2 y,xy 2 ,y 2 , 

x 3 y, x 2 y 2 , xy 3 , x 3 y 2 , x 2 y 3 , x 3 y 3 ], 
and (a(t)} i s the generalized bending coordinates 

j. 

{a} = [ai,a2i a 3 ) a 4,a5i a 6,a7,«8, a 9 5 a l0, 

ail,ai2>ai3,ai4,ai5,ai6] 

If the vector of the 16 bending nodal displacements is given as 


(B.2) 


(B.3) 


rp 

{u>b} = [wi,W x l,W y l,W xy l,W2,W x 2,lV y 2,Wx y 2, 


W3,W x 3,W y 3, W xy 3,W4, W x 4 , W y 4 , W xy ±) 


(B.4) 


then equation (B.l) and its derivatives at the four element nodes yield a set of 16 linear 
equations. 


{<«»} = pr i w 


(B.5) 
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which gives the generalized coordinates from the nodal bending displacements, so does 
the bending deflection distribution in the element, as 

M = [T b ]{w b } (B.6) 

where the transformation matrix between transverse displacements and generalized bend- 
ing coordinates, [T b ], is only a function of the coordinates of the four element nodes. 

The transformation matrix between inplane displacements and generalized inplane 
coordinates can be determined in a similar way. The inplane displacements are assumed 
to be bilinear, that is 

u(x,y,t) = b\ + 1>2X + b 3 y + 6 4 xy = [H u (x, y)]{6(f)} 

(B.7) 

v(x,y,t) = b 5 + 6 6 x -1- bjy + hxy = [H v (x, y)]{6(*)} 
where [H u (x,y)] and [ H v (x,y )] are the inplane shape functions 

[H u {x,y)\ - [l,x, y,xt/,0, 0,0,0] 

(B.8) 

[tf„(x,y)] = [0,0,0, 0, l,x,y,xy] 
and {b(t)} is the generalized inplane coordinates 

{b} T = [61,62, 63, h, h,be, 67, M (B.9) 

If the 8 inplane displacements at four element nodes are given as 

{u>m} = [wi,Ul,U2,y2,W3,U3,«4,U4] (B.10) 

Then equation (B.7) at the four element nodes yield a set of 8 linear equations, 

{ w m} = [Tin] {6} (B .11) 

which gives the generalized coordinates from the nodal inplane displacements, so does 
the inplane displacements distribution in the element, as 

{6} = [T m ]{w m } (B.12) 

where the transformation matrix between inplane displacements and generalized inplane 

coordinates, [T m \, is also only a function of the coordinates of the four element nodes. 
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B.2 Transformation between Nodal Displacements and Strains 

The transformation between the nodal displacements and strains can be determined 
directly from the strain — displacement relation. Substitute equations (B.l) and (B.7) into 
equation (2.22), we obtain the matrix in equation (A.4) as 

' - d 2 H w /dx 2 

[Cb]= - d 2 H w /dy 2 (B.13) 

.-2 d 2 H w /dxdy_ 
and 

'dH u /dx 

[Cm] = dH v /dy (B.14) 

pH u /dy + dH v /dx 
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APPENDIX C 

Coupling between FEM and BEM 


The velocity vector of each element in a constant rectangular BEM model is the 
displacement at the center of the element 

{u n } = [u on } T (C.l) 


while in FEM, the displacement vector is 

= [iCj, Wj x , W{y, Wixy, U,', Uj] 

T (C.2) 

{w} = 

If the mesh of FEM model and BEM model on the internal surface of the plate are 


identical, then from the displacement shape functions, we have. 


' H w (x 0 ,y 0 ) 

0 

fra | 

0 

H v (x 0 ,y o ) 

Ur m ]J 

0 

H v (x 0 ,y 0 ) 


(C.3) 


where 


x 0 = 
Vo = 


Xl + X2 + ^3 + £4 

4 

yi + n + yi + y± 

4 


(C.4) 


with xi and y, are the coordinates of element nodes. 

Due to the compatibility of the velocity at the inner surface of the plate, we have 
the element in vector {u n } as 


Une = {Pe} T [T e ]{«>} 


(C.5) 


and 

{&} = k nj,,n z ] T (C.6) 

where n x , n y , and n z are the components of the unit normal vector of the element plane. 
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